Entanglement Spectra and Entanglement 
Thermodynamics of Hofstadter Bilayers 



John Schliemann 

Institute for Theoretical Physics, University of Regensburg, D-93040 Regensburg, 
Germany 

E-mail: j ohn . schliemann@physik . uni-regensburg . de 

Abstract. We study Hofstadter bilayers, i.e. coupled hopping models on two- 
dimensional square lattices in a perpendicular magnetic field. Upon tracing out 
one of the layers, we find an explicit expression for the resulting entanglement 
spectrum in terms of the energy eigenvalues of the underlying monolayer system. 
For strongly coupled layers the entanglement Hamiltonian is proportional to the 
energetic Hamiltonian of the monolayer system. The proportionality factor, however, 
cannot be interpreted as the inverse thermodynamic temperature, but represents 
a phenomenological temperature scale. We derive an explicit relation between 
both temperature scales which is in close analogy to a standard result of classic 
thermodynamics. In the limit of vanishing temperature, thermodynamic quantities 
such as entropy and inner energy approach their ground-state values, but show a 
fractal structure as a function of magnetic flux. 



1. Introduction 

The study of quantum entanglement has by now developed to a mature subfield of 
many body physics [TJ [21 E]. Among the recent developments, the concept of the 
entanglement spectrum [I] has been applied to a plethora of different systems. These 
comprise quantum Hall monolayers at fractional filling [H El El El El QUI [HI H21 
[131 [TH [151 EH [TF1 [THl [19] where various spatial geometries and ways of separating 
the system into subsystems have been investigated. In the case of bilayer quantum 
Hall systems, a very natural way of defining a subsystem is provided by the double 
layer structure. Specifically, a numerical investigation of quantum Hall bilayers at 
total filling factor of unity revealed a striking similarity between the entanglement 
spectrum of the composite systems and the energy spectrum of a monolayer [20J. A 
analogous observation was reported slightly earlier in a numerical study of two-leg spin- 
1/2 ladders [2T] which was subsequently understood in terms of perturbation theory in 
the limit of strong rung coupling [221 [23], a result which is remarkably also valid for 
arbitrary spin length [21]. Similar findings were obtained in the simple case of coupled 
chains of free fermions where explicit expressions for the entanglement spectrum at 
arbitrary coupling can be derived [2"2l |2"o] . The starting point of the present work is 
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to extend these results to Hofstadter bilayers, i.e. coupled hopping models on two- 
dimensional square lattices in a perpendicular magnetic field [26J. The energy spectrum 
of such systems generates, as a function of the magnetic flux per unit cell, highly self- 
similar and visually appealing structures known as Hofstadter butterflies. As we shall 
see below, these features immediately translate to the entanglement spectrum via an 
explicit formula for the entanglement levels in terms of the energy eigenvalues of the 
monolayer system. Another focus of the present work are thermodynamic properties of 
the reduced density matrix and its entanglement Hamiltonian. In particular, we derive 
a thermodynamic entanglement temperature and inner energy starting from an effective 
coupling parameter of the bilayer system which can be viewed as a phenomenological 
entanglement scale. 

Entanglement spectra of Hofstadter monolayers were already investigated in 
Ref. [27] with the underlying square lattice being partitioned into two blocks. This 
and related ways of defining subsystems were also used in other studies of entanglement 
spectra of spin systems |2SlE3ll3nil^E2lE3ll31EaE5ll3Zll3SlESlSQlllI]^ Other recent 
investigations in connection with entanglement spectra include topological insulators 
[121 US], rotating Bose-Einstein condensates [44], coupled Tomonaga-Luttinger liquids 
[45] . interacting bosons [461 S3 HB] , and complex paired superfluids ]49j. 

This paper is organized as follows. In section [2] we review and extend previous 
results on free fermions in two-party lattice systems in the absence of a magnetic 
field. In the following section [3] we analyze the entanglement spectrum of Hofstadter 
bilayers and derive an explicit formula for the entanglement levels as a function of the 
underlying (highly fractal) monolayer system. Section[4]is devoted to a thorough analysis 
of the thermodynamic properties of the reduced density matrix and its entanglement 
Hamiltonian. We close with a summary and an outlook in section [5j 

2. Free Fermions without Magnetic Field 

Let us first briefly review and generalize results of Ref. [22] on free fermions on coupled 
lattices. We consider the Hamiltonian 



hopping term proportional to t± leading, at given wave vector, to a 2 x 2 eigenvalue 
problem whose elementary solution is 



U = J2 [£A{k)a±a % + e B $)bp % - t ± ( a +6g + b^) 



(1) 



k 



where a^,a^ (b$,b%) generate and annihilate fermions with wave vector k and energy 
EA{k) (£b(/c)) on a (i-dimensional lattice A (B). The two subsystems are coupled by a 




(2) 



k 



with 



e±{k) = \ (e A {k) + e B {k)) ± ^ At\ + (e A {k) - e B {k)) 



2 



(3) 
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and 



Pi 





t± 


4 a t - 


\t±\ 




t± 


\ a \ + 


\t±\ 



4 h t 



where 



\ 



1 ± 



e A (k) - e B (k) 



\ 



Let us now focus on the case where the above dispersion bands do not overlap, i.e. 



(4) 
(5) 

(6) 

(7) 
(8) 

holds identically in k since then we have £+(ki) > > e^{k 2 ) for all ki, k 2 . Indeed 
previous work [221 ES] has concentrated on the situation e A {k) = — s B (k) where the 
inequality ^ is obviously valid. For a half- filled system fulfilling ([7]), only the single- 
particle states generated by f3i are occupied in the ground state, and upon tracing out, 
e.g., subsystem B one obtains, following the very general arguments given in j5(]l EI], a 
reduced density matrix of the form 

Prcd = 77 exp(-H ent ) (9) 
with Z = tr(exp(— "H cnt )) and an entanglement Hamiltonian 



e+(h) >e-(k 2 ) 
for all wave vectors k%, k 2 . This is in particular the case if 
(e A (k)+e B (k)) 2 < (e A {k)-e B {k) 



(10) 



The entanglement levels £(fc) are given by 

In 1 V * 



(%)' 



In 



2 arsinh 



4tj + (e A (k) - e B (k)y + (e A (k) - e B (k)) 
At 2 ± + (e A (k) - e B (k)y - (s A (k) - e B {k)), 
e A (k) -e B (k) 



2\t ± \ 



(11) 
(12) 
(13) 



The exponential form ^ of the reduced density matrix can be derived by demanding 
that this expression should generate the same one-particle correlations in the monolayer 
subsystem as the underlying pure bilayer ground-state density operator. Using 
elementary properties of fermionic systems leads then to the explicit result (11) for 
the entanglement levels [50]. We note that the entanglement spectrum as discussed is 
generated by the single-particle operator ( 10 ) as appropriate for free fermions. This is 
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different from the situation in interacting systems where the eigenstates of the reduced 
density matrix are in general nontrivial many-body states. 

The above considerations extend the results of Peschel and Chung (22] to arbitrary 
dimension d > 1 and independent dispersions EA(k) ^ —esik) in the two subsystems. 
We note that the entanglement spectrum (13) depends, differently from the energy 
spectrum (jij), only on the difference (A;) — £#(&), but not on both quantities separately. 
Moreover, both the energy and the entanglement spectrum are invariant under a change 
of sign of t±. 

If the two energy bands overlap, i.e. condition ([7]) is violated, the situation cannot 
be analyzed in such general terms. Here the £(k) will in general not have support in 
certain parts of the first Brillouin zone, while in other parts two branches of entanglement 
levels will occur. 



3. Hofstadter Bilayers: Entanglement Spectra 



Let us now concentrate on two coupled two-dimensional square lattices with the 
Hamiltonian 



where 



Ha tA ^ ] ^m+l,n"' m ) n ®'m,n-\-l®" m i' n > h-C 

m,n 

W-B = - ts X] [bm+l, n bm,n + &m,n+l&m,n + h.C. 



and 



aib^ + h.c. 



~t± [ a m,nbm,n + h.C. 



(14) 

(15) 
(16) 

(17) 



Here a+ n , b^ nn generate particles at sites m,n in layer A and B, respectively, i.e. 
compared with equation ([T]) we have 

£A/B{h) = -2t A / B (cos(k x a) + cos(k y a)) (18) 

where a is the lattice constant in both layers. In order to implement a perpendicular 
magnetic field we work, as common, in Landau gauge with the discretized vector 
potential A = (0, Bma, 0) and concentrate on rational values of the magnetic flux per 
unit cell $ = Ba 2 /(h/e) = p/q (in units of h/e) where p and q are integers without 
any common divisor except unity. Thus, the system is periodic in the x-direction with 
periodicity qa while the periodicity in the ^-direction is still the lattice constant a. The 
standard Peierls substitution [52] leads the following ansatz for the state vector, 

1 



\k) 



N X Ny W.I! 



|^exp \ ik x ma + i yk y + —Bmaj na 
■ (a m (k)a+ jn + (3 m {k)b^ n )] |0) 



(19) 
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with a m (k) = a m+q (k) and j3 m (k) = (3 m+q (h). Here |0) is the fermionic vacuum, and 
N x , N y are the total numbers of unit cells in x- and ^/-direction, respectively, assuming 
periodic boundary conditions commensurate with q. The stationary Schrodinger 
equation 

H\k)=e{k)\k) (20) 
leads to a Harper equation being equivalent to the 2q x 2q eigenvalue problem 



H(k)v{k) = e(k)v(k) 
for the 2g-component spinor 



v(k) = (a (k), . . . , a q -i(k), /3 (k), . . . , P q -i(k) 



and 

H(k) 
Here the matrix 



h(k) 



t A h(k) -t±l qxq \ 
-t±l qxq t B h(k) J 

( r (k) z(k) 
z*(k) n{k) z(k) 
z*(k) r 2 (k) z(k) 

: 



(21) 
(22) 

(23) 



z*(k) \ 





z(k) 



V z{k) 







(24) 



z*(k) r q ^{k) J 



with r m (k) = 2cos(27r$m + k y a) and z(k) = exp(ik x a) corresponds to the classic 
Hofstadter monolayer problem [26J . The diagonalization of the latter quantity in general 
requires numerics, and explicit analytical results are possible only for very special values 
of the magnetic flux $ = p/q. 

Now let the unitary matrix u(k) diagonalize h(k), i.e. 



u(k)h(k)u + (k) = diag ^£0(^)5 ■ • • , £ q-i(k)) ■ 



(25) 



Since the off-diagonal elements in (23) are proportional to the q x q unit matrix, they 
remain unchanged upon a simultaneous diagonalization of both diagonal blocks. Thus, 
all four blocks are rendered diagonal, and the diagonalization of (23) is reduced to 2 x 2 
problems of the form 

/ t A e m (k) -t x 
\ -t ± t B e m (k) 
with the obvious eigenvalues 

1 



(26) 



\{k)=^-{t A + t B )e m {k)±^ltf 



\J±tl + ({t A -t B )e m {k)y 



and corresponding eigenvectors 



It^Vmik) 



(27) 



(28) 



\t±\ 
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Figure 1. Left panel: Classic monolayer Hofstadter spectrum (25) as a function of 
the magnetic flux. Right panel: Energy spectrum e+ (blue) and e~ (red) according 
to Eq. (271 of a bilayer system for t A = ts = t± . The spectrum consists of two copies 



of the Hofstadter butterfly shifted by ±t±. 



where 



Vm 



(*) 



\ 



it A ~ t B )e m (k) 



K \Jtt\ + ((t A - t B )e m (k)) 2 j 



(29) 



Fig. [T] shows the numerically computed classic Hofstadter spectrum of a single 
square lattice as a function of the magnetic flux $ per unit cell, along with the energy 
spectrum of a bilayer system for the particularly simple case t A = t B = t±. Here 
the Hamiltonian is invariant under exchange of layers, and the eigensystem consists of 
states being either symmetric or antisymmetric under this operation. Thus, as seen in 
the left panel of Fig. [TJ the bilayer energy spectrum comprises two monolayer Hofstadter 
butterflies shifted by ±tj_. 

We now concentrate again on the case where both groups of dispersion branches 
do not overlap, 



(30) 



for all fci, &2 and mi, m% G {0, . . . , q — 1}. Analogously to rt8J) the stronger condition 



(t A + t B f <(t A -t 



B 



(i.e. t A and t B need to differ in sign) implies e+ (&i) > > £ m . 2 (k 2 



(31) 

and therefore the 



inequality (30). Thus, under condition (30) only the single-particle states with energies 
e~ (k) are occupied in the ground state of a half- filled system. Tracing out again layer 
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B leads to entanglement levels of the form 



2' 



Uk) = In | 1 , | (32) 

ln | ' ySj + ((t A - t B )e m {k)Y + (t A - t B )e m {kj 
kt\ + ({t A - t B )e m (k)) 2 - {t A - t B )e m {k) J 



2a r sinh| ^-^ j. (34) 
These entanglement levels enter the entanglement Hamiltonian 

9-1 

^ent = J2 E^(^) a L a fc™ ( 35 ) 
m=0 £ 

at given flux $ = p/q. Here the operators ai" m generate Hofstadter monolayer states 
with eigenvalue e m (k), i.e. the monolayer Hamiltonian including the magnetic field can 
be formulated as 

9-1 

n A = t A M*04m a fc™ • ( 36 ) 



m=0 



k 



We note that the matrix ( 24 ) has the obvious property 

h(k) = -h(k + tt) (37) 
with 7? = (it /a, ir /a) which implies 

Em(fc) = - e,-l-m(* + 7?) > (38) 

e±(fc)= - e T_ 1 _ m (fc + ^), (39) 

£m(£) = -^-l-m(fc + 7T), (40) 

where we have the eigenvalues e m (k) of /i(fc) assumed to be given in ascending order. 
These relations enable to use, with some caveat and specifications, the expression (34) 
for the entanglement levels even in situations where the condition (30) is not fulfilled: 
The left panel of Fig. [i] shows the dispersions e m (k) for a magnetic flux of $ = 1/7 
and typical parameters. As seen in the figure, the eigenvalues form rather flat bands 
which do not overlap. The right panel displays the corresponding bilayer energy s^ik) 
for parameters violating the inequalities (30), (31). In the ground state of a half-filled 
system the lowest q = 7 bands are completely occupied such that the highest occupied 
(valence) band is ef(k), which is also the highest band among the branches e^(fe) 
having negative energy. Note that ef(k) does also not overlap with any other band 
e^(fc); otherwise such two bands would only partially be filled each. Let us therefore 
focus on this "insulator" scenario where the valence band of a half filled system does 
not overlap with other bands implying that in the ground state all bands are either fully 
occupied or empty. Let us now assume that some band e+(fe) has negative energy and 
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Figure 2. Left panel: Eigenvalue dispersions £ m (fc) of the Hofstadter monolayer 
for $ = 1/7. Right panel: Corresponding bilayer energy dispe rsio ns £ + (fc) (blue) and 
e~ (fc) (red) for a choice of parameters violating the conditions (30 1, (31 ). The situation 
for other fluxes $ = p/q is qualitatively similar, but with increasing q more difficult to 
display due to the larger number of bands. 



is therefore fully filled. Thus, according to Eq. (39) £„-i- m (k) has positive energy and 



is therefore empty. The entanglement branch arising from e^k) is 

In 



-1— m 



+ 7?) 



(41) 



where we have used Eqs. (29),fl39h. Thus, up to a rigid shift in the wave vector argument, 



the entanglement levels arising from the occupied band s^JJk) reproduce the missing 
branch corresponding to eg_!_ m (fc). In this sense (i.e. with the wave vector dependence 



being suppressed) the expression (34) for the entanglement spectrum can also be used 



even if the inequality (30) does not hold, but all bands are, in the ground state, either 



empty or completely occupied. On the other hand, starting directly from the original 



2q x 2q eigenvalue problem (21), the entanglement levels can be expressed as 



im{k) = In 



Cm(k) 



Cm(k) 



where 



9-1 



Cm(k) — ^ 



\a 



(m) 



(fc) 



(42) 



(43) 



m'=0 



and the a^\k) are the components referring to subsystem A in the eigenvectors (22) 



corresponding to the lowest q eigenvalues in (21 ), m e {0, . . . , q — 1}. Where applicable 



see above), the results (34) and (42) of course coincide, as we have also checked by 



explicit numerical calculations. 
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Figure 3. Energy spectrum e+ (blue) and e m (red) according to Eq. (27 1 of a bilayer 
system for tj± = ij_ and various values of t/3. 



In summary, the expressions (27) and (34) provide, in full analogy with 



Eqs. ([3J),( 13 ), explicit relations between the energy and the entanglement spectrum, 
respectively, of the composite system and the energy spectrum of the monolayer. 
Moreover, the entanglement spectrum £ m (k) depends, for a given monolayer spectrum 
£ m (k) only on the single quantity 
tA — ts 



A 



\t±\ 



(44) 



which will in the following be referred to as the (effective) coupling parameter. 

As already discussed above, in the case A = the layers are maximally entangled 
with each other, and the all entanglement levels are zero. In Fig. [3] we have plotted 
the energy spectrum of the bilayer system according to Eq. ( J27| ) for tA = t± and 
various values of ts- As seen there, if £4 and ts differ in sign (here: ts < 0), an 
energy gap which is independent of the magnetic flux $ occurs. Fig. [4] displays the 
corresponding entanglement spectra which form Hofstadter butterflies "deformed" by 



the inverse hyperbolic sine function occuring in Eq. (34). 
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Figure 4. Entanglement levels of bilayer systems for the same parameters as in 
Fig. [3] The spectra are Hofstadter butterflies "deformed" by the inverse hyperbolic 
sine function occuring in Eq. (34). The effective coupling parameter A is defined in 
Eq. d44|. 



4. Hofstadter Bilayers: Entanglement Thermodynamics 

In the spirit of Ref. [20] we now investigate the thermodynamics based on the 



entanglement Hamiltonian (|35j) and the pertaining reduced density operator p re d- From 

(45) 



Eq. (34) one finds 

U(k) = \e m (k) - y (\e m (k)) 3 + O (A 5 ) 



implying (cf. Eq. (36))) 



n 



ent 



(46) 



for small coupling parameter A. This suggests to interpret A as an inverse "entanglement 
temperature" and approximate the reduced density matrix as p ro d ~ Po with 

Po = v exp(-\n A /t A ) (47) 
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Figure 5. The entropy s(A, $) per unit cell as a function of the magnetic flux $ for 
various values of A. For A 10 and larger, fractal features of s(A, $) become visible. 
The data (as well as the one in the following Figs. [6] and [7]) was calculated numerically 
according to Eqs. (51 1,(52 ) using an appropriate discretization of the magnetic Brillouin 



with Z = tr(exp(— A"HaAa)) and "Ha/^a describing energy of the subsystem. The 
above findings are completely analogous to the ones in Ref. [20] on quantum Hall bilayers 
where the above approximate relation (46), valid for strongly coupled layers, was used 
to numerically analyze a quantum phase transition in the total double layer system. 
However, having a closed analytical result for the entanglement levels at hand, we shall 



follow a different route here and employ the full expressions (34),(j35j) (not necessarily 
approximated by their first order in A) to evaluate thermodynamics in a more systematic 
manner. As we shall explore in more detail below, the coupling parameter A can be 
used as a phenomenological temperature scale, but it is not identical to the inverse 
thermodynamic temperature given by the derivative of the entropy with respect to the 
(appropriately defined) inner energy. 

A very simple issue is the average particle number iV = (N) with iV = J2 m k a \ m a km 
and (■) = tr(p roc r). Here one always has, independently of coupling parameter and 
magnetic flux, iV(A.$) = N x N y /2 which is clear from the fact that the total bilayer 
system is half-filled and unbiased with respect to its subsystems. Formally this result 
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Figure 6. The entropy s(A, $) per unit cell as a function of the magnetic flux <f> for 
A ^> 1. The figure contains all fluxes $ = p/q with p,q coprime and q < 100. 



can be established, in the thermodynamic limit N x ,N y 3> 1, as follows, 



iv(A,$) = i\yv 3 



J/ (2vr) 2 



9-1 

E 

m=0 



27r/(qa) 



2-ir/a 



elk*. 



I _|_ g€m(fc) 



N x N y - 
N x N y /2, 



9-1 
1 m=0 



27r/(ga) 



2ir/a 



i 



+ 



Ae m (fc)/4 



l + (Ae m (A;)/2) 2 , 
(4£ 



where we have used Eqs. (33) and (38). This constancy of the averaged particle number 
can of course also be seen as a consequence of the fact that a reduced density matrix 
of the form ^ can be viewed as a grand-canonical statistical operator with constant 
chemical potential fi = 0, which lies exactly in the middle of the symmetric entanglement 
spectrum. 

Let us next turn to the other extensive quantities entropy and energy which we 
define, again in the thermodynamic limit, by 

S — (— lnp red ) = N x N y s(X, $) , (49) 
E = (H ent ) = N x N v e(\,$) (50) 
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Figure 7. The energy e(A, $) per unit cell defined by Eq. (52) as a function of the 
magnetic flux $ for A ^> f . 




Figure 8. Entropy s(A, $) (left panel) and energy e(A, $) (right panel) per unit cell 
as a function of A for various fluxes The inset shows the derivative ds/de. All 
quantities decrease monotonously with increasing A and do not show any apparent 
qualitative difference for different $. 
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with 



2 9-1 ,27r/( 9 a) r 2Tv/a 111 f 1 + e U ^ 

AU I AU V 



S ( A > $ ) = 7^^E/ dk * dk y^ T7^' ( 51 ) 

(27r) i j ^ Jo Jo \ + e U(k) 

„2 9-1 r 2Tv/{qa) r 2-K/a £ (t\ 

e(A, d>) = £ / dfc, / ^rff^ • (52) 

(53) 

The bar at the above quantity E is meant to indicate that the definition of energy will 
be subject to some refinement further below. Moreover, as a consequence of Eq. (40), 
this quantity is non-positive, 

e(A,$)<0 (54) 

for all A and $. Besides, since the entropy S is obviously proportional to the area of 
the system coinciding with the boundary to the other subsystem, the so-called area law 
is fulfilled [3]. 

For small A 1 one obtains the expansions 

s(A, $) = \ In 2 - ^ + A 4 ^ (7 + 2 cos(2tt$)) + O (A 6 ) , (55) 

e(A, $) = - A 2 + ^ (7 + 2 cos(2tt$)) + O (A 6 ) , (56) 
where we have used the fact that 

„2 9-1 f 2n/(qa) r2ir / a , _ . n 

: = 7^ E / ^ / <fty (*(*)) (57) 

a 2 r2ir/(qa) r2n/a 

= (2^io ^io ^ y 
vanishes for odd n, and for low even values one has ujq(§) = 1, W2($) = 4, and 



(h(k)) n ] (5f 



W4 ($) = 4(7 + 2 cos(2tt$)) , (59) 



as it is readily derived from the explicit form (24) of the matrix h(k). We note that the 
zeroth-order value s(0, $) = (ln2)/2 of the entropy per unit cell is due to the fact that, 
at strong coupling, each particle ends up in either layer with equal probability, and the 
total bilayer system is half-filled. 

In Fig. [5] we have plotted the entropy s(A, $) per unit cell as a function of the 
magnetic flux $ for various coupling strength starting with A being of order unity, 
where s(A, $) is a smooth function of $. At values of A ~ 10 and larger, however, 
s(A, $) becomes more irregular and shows typical fractal features such as self-similarity. 
As seen in Fig. [6j this property of s(A, $) becomes more pronounced with increasing A. 
A similar behavior is found for e(A, $), as illustrated in Fig. [7j 

Fig. [8] shows the entropy s(A, $) and the energy e(A, $) as a function of A for a few 
representative values of $. As seen, both quantities decrease monotonically with A, and 
in this sense A qualifies as a phenomenological (inverse) temperature scale. However, if 
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A were the true inverse thermodynamic temperature, standard thermodynamics would 
require it to equal the derivative 

ds 0s(A,$) /de(A,$)\ _1 



de OX \ OX J 
where the derivatives with respect to A on the r.h.s can be expressed as 

<9s(A,$) a 2 ^ rWM r**fa 

/ 4 I dk x I dky 



(60) 



d\ (2tt) 2 



m=0 



(- In (l + e 6 *®) + l) 



2 



(61) 



(62) 



(l + e M*)) (l + e -U(*)) y/l + (Xi(k)/2) 

de(X, $) _ a 2 ^ /-2t/(9«) 

<9A _ (27T) 2 ,~ To dky 

(-em(fc) + (l + e^))e(fc) 

(l + e M*)) (l + e-^W) y/l + (\e(k)/2) 

In particular, for small A 1 one finds from Eqs. (55) and (56) 

independently of the magnetic flux $. As seen in the inset of Fig. |8j and also very 
explicitly by Eq. (63), the derivative (60) is certainly not equal to A. 

The reason for this behavior is that the reduced density matrix formulated as 
Pred = exp(— l-L eiL t}/Z does not match a canonical equilibrium state characterized by an 
inverse temperature /3. Let us therefore rewrite the entanglement Hamiltonian according 
to 

?4 nt (A, $) = /3(A, $)-H C an(A, $ ) , (64) 

where the inverse thermodynamic temperature /3(A, $) is determined as a function of A 
as follows: Defining the thermodynamic inner energy 

£(A,$) = (H can (A,$)} (65) 

and the free energy 

F(A,$) = £(A,$) - 5(A,$)//3(A,$) (66) 
= -ln(Z(A,$))//3(A,$) (67) 

one easily verifies 
dF 

"w = E m 

with F = (3F = E — S, where we have used the stipulated relation 

(69) 
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The above equations (68) and ( p9J ) are indeed equivalent, and from Eq. (68) it follows 

d\n/3 1 dF _ ld(e-s) 
E ~d\ ~ i 



(70) 



which is the sought equation determining /3(A). For small A the expansions (55) and 

(71) 



(56) lead to 



din /3 



A + 35 A (7 + 2cos(27r$)) + 0(A 2 ) 



such that 



m/3(A,$) 



In k + - In A + 



(72) 



2 ... 24 A 2 (7 + 2cos(27r$)) + 0(A 3 ) . 

Here the integration constant Ink, k > 0, reflects the freedom of choosing a unit 
to measure /3, i.e. k plays the same role as Boltzmann's constant in standard 
thermodynamics. The entropy s and energy e = e//3 per unit cell read as a function of 
small 



1 



In 2 



4 



(f3/k) 4/3 + O (/3 3 ) 



e(/3,$)= -i(/3A) 1/3 + G(/3 2 



(73) 
(74) 



Remarkably, the inverse thermodynamic temperature 

/3(A) = A;A 3/2 + O (A 7/2 ) (75) 
scales for A 1 as A 3 ' 2 and is not linear in this parameter, differently what would follow 



from the ansatz p re d = p , cf. Eq. (47). The technical reason for this observation is that 



Po is, while being of canonical form, not the linear approximation to p rc d but contains 
also arbitrary high powers in A. 

At large A, i.e. weak coupling between the layers, the entropy s characterizing their 
mutual entanglement will (along with its derivative) eventually vanish. Thus, for A ^ 1 



Eq. (70) can be simplified as 
din (3 d In lei 



leading to 



/3(A,$) w fc|e(A,$)| = -fce(A,$) 



(76) 



(77) 



where we have used Eq. (54) and have adjusted the constant k consistent with Eq. (72). 



In particular, the above relation implies 

1 



e(A,$) 



(7E 



for large A. We note that the sign on the r.h.s. of Eqs. (78) depends on the inequality 



(54) and would be different for positive e(A, $) at A ^> 1. 

The inset of Fig. [9] shows (3 as a function of A obtained via numerical integration 



of Eq. (70) for the same flux values $ as in Fig. 8 The data depends only weakly 



on $ and follows the power law (75) at small A while the behavior at large A is well 
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Figure 9. Entropy s (left panel) and energy e (right panel) per unit cell as a function of 
the inverse thermodynamic temperature /3 for the same flux values $ as in Fig. [8j The 
inset shows j3 as a function of A. This data was obtained via a numerical integration 



of Eq. ( 70 ) starting from expression ( 72 1 with an appropriately small but finite initial 



argument < A <C 1. 



described by a logarithmic dependence. In the main panels of Fig. [9] we have plotted 
the entropy s(/3, $) and the energy e(/3, $) (as opposed to e) as a function of (3. for 
large (3, e(/3, $) levels off and converges (slowly) to e = — 1/k, consistent with Eq. (78). 
Such a saturation of the thermodynamic energy should be expected since both layers 
get more and more decoupled in this limit. The originally defined energy e = /3e does, 
in accordance with Eq. (77), not show such a behavior but decreases unboundedly with 
increasing A, as shown already in Figs. [7] and [8] While approaching their ground state 
values in the limit A — > 00 the inverse thermodynamic temperature /3(A, <&) as well as 



the inner energy e(A, $) per unit cell display fractal features, as illustrated in Fig. 10 
An analogous behavior was found for the entropy in Fig. [6j 

In summary, the parameter A can be viewed as a phenomenological inverse 
temperature scale having intuitive properties like the decrease of energy and entropy 
with increasing A. The thermodynamic inverse temperature fulfilling standard 
thermodynamic relations, however, is given by /3(A). Eq. (70) establishes a 1-to-l 
mapping between both quantities and is completely universal in the sense that its 
derivation does not depend on any detail of the underlying system or its entanglement 
Hamiltonian. Indeed, a relation very similar to Eq. (70) (with in fact an identical 
l.h.s.) is obtained in standard thermodynamics when connecting the thermodynamic 
temperature to a phenomenological temperature scale obtained from, say, the change 
in volume of a given body upon changing its internal energy [53j El] • Also the results 
(76) to (78) are (up to the sign (54) of e) very general since they only rely on the 
vanishing of the derivative of the entanglement entropy s(A) in the limit of weak coupling 
A ^> 1 /3 ^> 1). This statement is of course just the analog of the third law of 
thermodynamics which requires the entropy to approach a constant in the limit of zero 
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Figure 10. The inverse thermodynamic temperature /3(A, $) (upper two panels) and 
the inner energy e(A, <£>) per unit cell (bottom panel) as a function of the magnetic 
flux $ for large A ^> 1. The figure contains all fluxes $ = p/q with p,q coprime and 
q < 100. 



temperature. In the same limit, the thermodynamic energy e per unit cell reaches 
according to Eq. (78) a saturation corresponding to the ground state energy in classic 
thermodynamics. Remarkably, within the formalism of entanglement thermodynamics 
outlined here, this limit is universal and depends only on the constant k, i.e. the unit 
chosen to measure temperature. 



5. Conclusions and Outlook 



We have derived an explicit expression for the entanglement levels of Hofstadter 
bilayers in terms of the energy eigenvalues of the underlying monolayer system. 
For strongly coupled layers we find the (expected) proportionality between the 
entanglement Hamiltonian and the energetic Hamiltonian of the monolayer system with 
the proportionality factor given by an effective coupling parameter. This parameter, 
however, is not identical to the inverse thermodynamic temperature, but represents a 
phenomenological temperature scale. We have devised an explicit relation between 
both temperature scales which is in close analogy to a standard result of classic 
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thermodynamics. The introduction of the thermodynamic temperature also implies a 
redefinition of the inner energy. In the limit of vanishing temperature, thermodynamic 
quantities such as entropy and inner energy approach their ground-state values, but 
show a fractal structure as a function of magnetic flux. 

The relation between the phenomenological temperature scale (given by an 
appropriate coupling parameter of the total system) and the thermodynamic 
entanglement temperature applies certainly also to other and more general situations. 
For instance, in Ref. [20] an entanglement temperature was obtained for bilayer quantum 
Hall systems by a numerical fit to exact-diagonalization data as a function of layer 
separation. In the light of the present work, this temperature scale should be seen as a 
phenomenological one. The analysis of the quantum phase transition performed there, 
however, should qualitatively not be affected by this issue, since the relation between 
a phenomenological and the thermodynamic temperature is expected to be a smooth 
function. On the other and, very explicit investigations as done here for non-interacting 
particles are obviously more difficult for an interacting system since clearly less analytical 
tools are available [20J. 

Similar comments apply to the situation of spin ladders in the limit of strong rung 
coupling J22J [23J, El]: The proportionality factor found there between entanglement 
spectrum and energy spectrum should also be viewed as a phenomenological temperature 
scale, but not the thermodynamic one. Moreover, Ref. [55] also introduces, studying 
block entanglement in a spin-1/2 XX-chain, an effective inverse temperature which 
grows monotonically with increasing block size. Although the situation there and the 
present one show obvious differences their possible interrelations are of interest. 

Finally, another possible extension of the present study is to consider other types of 
lattices. Given the large deal of interest devoted presently to graphene, the hexagonal 
geometry [SB] and its bilayer versions [57] are obvious candidates . Indeed, it is an 
interesting speculation whether, for example, the exponent 3/2 occuring in Eq. (75) 
depends on the lattice geometry. 
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